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BOUNDARY- LAYER TRANSITION AND DISPLACEMENT-THICKNESS 
EFFECTS ON ZERO- LIFT DRAG OF A SERIES OF 
POWER- LAW BODIES AT MACH 6 

By George C, Ashby, Jr., and Julius E, Harris 
Langley Research Center 

SUMMARY 

The influence of viscous effects on the zero-lift drag of a series of power-law bodies 
has been examined analytically and experimentally at a Mach number of 6. Bodies con- 
strained to length-diameter ratios of 6.63 with values of the power-law exponent in the 
shape equation of 0.25, 0.50, 0.60, 0.667, 0.75, and 1.00 were studied. The pressure drag 
of each body was computed by using the method of characteristics, and the skin friction 
was computed by using a numerical solution of the compressible boundary layer with 
the local flow properties predicted by the method of characteristics. The location of 
boundary-layer transition on the models was determined from phase- change -coating heat- 
transfer measurements, and the effects of transition and boundary-layer displacement 
were included in the calculations. The calculated drags were compared with experimental 
results obtained at a Mach number of 6 over a Reynolds number range, based on body 
length, of 1.50 x 10® to 9.50 x 10®. The analysis showed that the exponent for minimum 
drag with either inviscid or fully laminar flow was approximately 0.667. However, if 
boundary- layer transition existed on the bodies, the exponent for minimum drag was lower. 
For the experimental results at the higher Reynolds numbers, boundary-layer transition 
occurred, and the resulting exponent for minimum drag was approximately 0.6. 

The effect of displacement thickness on total drag was significant but was nearly 
constant for all the bodies exerting only a minor influence in the determination of the body 
with minimum zero-lift drag. For the present tests, the boundary-layer displacement- 
thickness effect on total drag was as high as 9 percent with the major contribution coming 
from wave drag. 


INTRODUCTION 

Early studies, typified by references 1 and 2, which were conducted to determine the 
body and/or wing profiles with minimum zero-lift drag considered only the inviscid drag 
and used approximate pressure laws. Subsequent investigators included the effects of 



skin friction but in a simplified manner, (See, for example, ref. 3.) The results from 
such basic investigations, however, have been useful design guides and have been shown 
to be generally applicable under lifting as well as zero-lift conditions (refs. 4 and 5). 

More recently, with the focus of attention on the advanced space transportation sys- 
tems, the practical application of body shaping to improve hypersonic performance has 
received greater impetus, and some refined results have been obtained. For example, in 
reference 6 a critical assessment of the pressure laws used in previous optimization 
studies (primarily the Newton pressure laws) revealed that more exact pressure laws, 
such as the method of characteristics, defined optimum bodies which had more volume 
forward than previous studies had shown. Furthermore, experimental data for a series of 
power-law bodies obtained at a Mach number of 6 over a large Reynolds number range in 
support of reference 6 indicated that a blunter body than the one predicted by the exact 
pressure laws would have the minimum zero-lift drag. This result was attributed to the 
variation of boundary- layer transition between the bodies at the test flow conditions and 
promoted the current study. 

The purpose of the present investigation was to determine the influence of viscous 
effects on the minimum zero-lift drag. Accordingly, calculations were made of skin- 
friction and boundary-layer displacement effects on a series of power-law bodies which 
were tested at a Mach number of 6 and reported in references 6 and 7. In the reference 
tests the condition of the boundary layer was not established, but the range of Reynolds 
numbers was large enough for the boundary layer to vary from laminar over the entire 
body to a laminar-transitional- turbulent state over the entire body. For the present study, 
the boundary-layer transition location was determined on reproductions of the force models 
by phase -change -coating heat-transfer measurements at zero angle of attack for Reynolds 
numbers corresponding to those of the reference test. (See refs. 6 and 7.) These meas- 
ured transition locations along with a pressure distribution calculated by the inverse 
method for the nose and the method of characteristics for the body were used in a 
nonsimilar-boundary-layer program to compute the skin-friction and wave-drag coeffi- 
cients on each body. An iterative approach was used to couple the boundary -layer and the 
pressure-distribution programs to account for boundary -layer displacement effects. The 
computed results are presented and compared with experimental drag coefficients herein. 

SYMBOLS 

Measurements and calculations were made in the U.S. Customary Units. They are 
presented herein in the International System of Units (SI) with the equivalent values given 
in the U.S. Customary Units. 
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A 


model base area, meters^ (feet^) 


drag coefficient, 


Drag 


total skin-friction-drag coefficient, 


Total skin-friction drag 


total-drag coefficient, sum of wave and skin-friction drag. 


Total drag 


wave -drag coefficient. 


Wave drag 


local skin-friction coefficient 


... . , Pe - 

pressure coefficient, ^ 


model length, centimeters (feet) 


exponent in shape equation for power-law bodies, ^ (see fig. 1) 

body surface pressure, newtons/meter^ (pounds/foot^) 
free- stream static pressure, newtons/meter^ (pounds/foot^) 
free-stream dynamic pressure, newtons/ meter^ (pounds/foot'^) 
base radius of power -law bodies, centimeters (feet) 


Reynolds number based on free-stream properties and model length 

coordinates on meridian curve of power-law body, centimeters (feet) 
(see fig. 1) 

distance along body surface, centimeters (feet) 

velocity at edge of boundary layer, centimeters/ second (feet/ second) 
end of boundary-layer transition 
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. beginning of boundary-layer transition 

6 inclination of body surface relative to free stream, degrees 

5* boundary-layer displacement thickness, centimeters (feet) 

MEASURED LOCATIONS OF BOUNDARY -LAYER TRANSITION 

Models 

The heat-transfer models used in this investigation were cast reproductions of the 
force models (ref. 7) and are shown in figure 1. The six power-law bodies, which are 
defined by ^ with n = 0.25, 0.50, 0.60, 0.667, 0.75, and 1.00, were circular in 

cross section, had a fineness ratio of 6.63, and were constrained for constant length and 
base diameter. Each model had an aluminum nose and a high-temperature epoxy plastic 
body. The measured coordinates of the heat-transfer models are within 0.040 cm of those 
for the force tests in reference 6. The noses of the heating models were individually 
shaped to match those of the force models. 

Wind Tunnel and Tests 

The transition measurements were conducted in the same tunnel as the force tests 
of reference 7, the Langley 20-inch Mach 6 tunnel. It is a blowdown type exhausting to 
either atmosphere or vacuum with an operating range of stagnation pressures from 3 to 
35 atmospheres (0.304 to 3.55 MN/m^) and stagnation temperatures up to 589 K. The 
general details of the tunnel, along with schematic drawings, are presented in reference 8. 

The location of boundary-layer transition was measured on each body, by using a 
phase -change -coating technique described in reference 9, at five Reynolds numbers, from 
1.50 X 10® to 9.50 X 10® based on model length. Four of the Reynolds numbers matched 
those of the reference force tests. Since transition did not occur on any of the bodies at 
the lowest Reynolds number, an additional intermediate Reynolds number was selected to 
give a more complete trend. All tests were conducted at zero angle of attack and sideslip. 

Method of Determining Transition 

Transition could be observed visually on these models because of heat-rate varia- 
tions between laminar, transitional, and turbulent boundary -layer regions. (See fig. 2.) 
Note that the light region in the photograph is the unmelted area. In the phase-change 
technique the heating rate is inversely proportional to the time required to melt the coat- 
ing, Progressing rearward from the apex of the model nose, the heating rate decreases 
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through the laminar region and becomes a minimum at the start of transition. The rate 
increases through the transition region until fully turbulent flow is established and then 
decreases at a relatively low rate in the turbulent region. Under these conditions, a coat- 
ing which melts at a specific temperature when applied uniformly to the model will begin 
to melt at the nose and then at or near the rear, depending on whether or not transition is 
completed on the body. The two melting regions approach each other and meet at the ini- 
tial transition point, which is the last place the coating melts. The series of photographs 
in figure 2 show an example of the melting sequence. 

For each test, the model was coated and mounted on a tunnel injection strut in the 
retracted position. After starting the tunnel and setting the flow conditions to the desired 
Reynolds numbers, the model was injected and observed visually until the phase-change 
process was complete. The time histories of the phase change of the coating were 
recorded with motion-picture cameras on both the top and right side of the models. The 
location of the beginning of transition (last melt point) and, where possible, the end of 
transition (first rearward melt point) were scaled from those pictures. 

Locations of Boundary- Layer Transition 

Typical peripheral distribution of the beginning of transition obtained from the phase- 
change- coating tests at o! = 0° are presented in figure 3, and the values used in the 
boundary-layer calculations are shown in figure 4. The average values were used because 
transition occurred at a smaller value of x. . on top than on the side and, in general, was 
not uniform about the body. Similar results were reported in references 10 and 11. For 
the present paper, the variation was less than 20 percent of the model length and could 
have been caused by any of several factors. For example, model irregularities could alter 
the location of boundary- layer transition. In reference 11 a deliberate nose-tip offset of 
0.1 percent of the tip length (0.005 cm) caused transition location to vary circumferentially 
about the body by 10 to 35 percent of the length. Nonuniform distribution of the phase- 
change coating or small flow angularity could also cause an asymmetric distribution of 
transition location. 

The variations in the location of the beginning of transition with Reynolds number 
for the five bodies in figure 4 are reasonably well behaved and do not exhibit any large 
anomalies. The apparent reversal of transition in the data for the bodies with n = 0.60 
and 0.75 at the higher Reynolds number is believed to be associated with nonuntformity 
of transition around the body. No reason is known for such transition reversal under these 
flow variations; therefore, the curve presented in figure 4 for n = 0.60 was faired with- 
out reversal. Calculations showed that, for the body with n = 0.60 at a Reynolds number 
of 9.5 X 10® and a transition- location variation from 52 to 69 percent of the body length, 
the skin-friction drag coefficient would vary only about 0.001, which is within the accuracy 
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of the measured drag coefficients. Thus, the difference between the curve as shown and 
a fairing through all the points was found to have an insignificant effect on skin-friction- 
drag coefficients. 

The distance to transition from the nose varies inversely with the power- law expo- 
nent (directly with nose bluntness, i.e., bluntness delays transition), which is in agreement 
with the results presented in reference 12. As the Reynolds number increases, transition 
moves forward, but the rate of forward movement is a direct function of n (inverse func- 
tion of nose bluntness}. At the higher Reynolds numbers of these tests, the rate of for- 
ward movement approaches zero for the bodies with n = 0.60 to 0.75; whereas, the loca- 
tion of transition on the body with n = 1.00 is more linear with Reynolds number. The 
difference in the variation of transition location between the body with n = 1.00 and the 
others at the high Reynolds numbers is not understood, but similar results were reported 
for a slightly blunted 5® half-angle cone at a Mach number of 8 in reference 13. 

The end of transition could not be positively identified for all cases (i.e., the begin- 
ning of turbulent flow did not occur on the body); those for which it could be identified are 
shown in figure 5. The ratio of the location of the end of transition to the location of the 
beginning of transition varies from approximately 1.2 to 1.7. Because no systematic vari- 
ation exists, a nominal value of 1.5, represented by the dashed line, was used for all cal- 
culations. A value of 1.5 for the average extent of boundary-layer transition is at the 
lower end of the range of values mentioned in reference 14 for a large amount of data. 

VISCOUS EFFECTS ON POWER-LAW BODY FOR MINIMUM 

ZERO- LI FT DRAG 

Effect of Boundary- Layer Transition 

The variation of the calculated skin-friction, wave, and total-drag coefficients with 
power-law exponent for the bodies at the four test Reynolds numbers are presented in fig- 
ure 6. The procedure used to compute the total-drag coefficients for the power-law 
bodies, including the effects of boundary-layer transition and displacement thickness, is 
presented in the appendix. Note that the measured locations of transition were used in the 
calculations. By considering wave drag only, the body with minimum drag would have an 
exponent n of approximately 0.667 (middle group of curves in fig. 6). This body is 
blunter than that predicted by Eggers et al. and by Miele (n = 0.75) in references 1 and 3, 
respectively, by using Newton’s pressure law ^Cp oc sin26). The difference in the bodies 
for minimum drag between this investigation and those of Eggers and Miele is attributed 
to the difference in pressure laws. This investigation used the method of characteristics, 
and Eggers and Miele used Newton's pressure law. A comprehensive discussion of the 


6 



effects of several pressure laws and methods used to predict minimum wave-drag shapes 
is contained in reference 6. 

When viscous effects are included ®)> exponent for minimum drag 

remains near 0.667 until the Reynolds number increases beyond 4.37 x 10 , above which 
it shifts towards a value of 0.60. This trend is caused by the variation of boundary-layer 
transition with Reynolds number for each power-law body (fig. 4). The location of transi- 
tion onset from the body nose ^ at a given Reynolds number varies inversely with the 
power-law exponent n (directly with bluntness, i.e., bluntness delays transition), while 
the rate of forward movement of transition location with Reynolds number dxt^i/dR«„^ 
varies directly with n. These divergent trends cause the increase in skin friction to be 
larger for the higher exponent bodies as Reynolds number increases and therewith the 
minimum total drag to shift towards the lower exponent, blunter bodies. Although these 
calculations include displacement thickness, it has little effect on the power-law body for 
minimum drag as is shown subsequently. 

The results of figures 4 and 6 confirm and explain more completely the shift in 
minimum drag towards a blunter body exhibited by the experimental data for the same 
series of bodies at a Mach number of 6 and the higher Reynolds numbers (refs. 6 and 7). 

In the references the shift was correctly attributed to the variation of the location 
boundary-layer transition from body to body at a constant Reynolds number. 

A comparison of the calculated values of total-drag coefficient with the experimental 
data is presented in figure 7. The height of symbols in the figure is representative of the 
experimental measuring accuracy, and the agreement between the measured and computed 
data for the majority of the points is within this accuracy. Of those points for which the 
measured and computed data do not agree, all but two are for the sharper bodies at the 
higher Reynolds numbers. H th.e location of transition for the bodies with n - 0.667 
and 0.75 had been the linearly extrapolated value (21 percent of length in fig. 4) and for 
the body with n = 1.00 had been at 50 percent of the body length, then the calculated val- 
ues of Cd T ^= 9 . 5 x 106 would have been represented by the dashed-line 

curves in figure 7. Under these assumptions the agreement between the calculated values 
and the measured data is good. At R^^^ = 6.42 x lO^, the calculated-drag- coefficient 
value for the body with n = 1.00 can be adjusted to agree with the measured value by 
assuming that the location and the nondimensionalized extent of transition are 58 per- 
cent of body length and 2.0, respectively, instead of the measured values of 46 percent 
and 1.5. Similar assumptions would bring the data for the body with n = 0.75 at the 
same Reynolds number into agreement with experiment. This exercise indicates that 
the calculated data agree with the experimental data for reasonable assumptions of 
boundary-layer transition. 
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Effect of Boundary-Layer Displacement Thickness 

Boundary-layer displacement thicknesses for each body are presented in figure 8 as 
a function of S/Z at each of the test Reynolds numbers. The location of transition used 
in the calculations is indicated on each curve if transition occurred for that particular 
Reynolds number. When transition begins somewhat forward of the body base, a marked 
dip in a* is followed by a rapid increase. However, the maximum thickness never 
exceeds that for the completely laminar-flow case at the lowest Reynolds number. The 
scale at the right (fig. 8) is a measure of the displacement thickness at the base in percent 
of the base radius. For the conditions of this investigation, the displacement thickness at 
the base varies from approximately 4 to 12 percent of the base radius. 

Effect on surface pressure . - The computed pressure distributions on each body with 
and without boundary-layer displacement effects at the highest and lowest Reynolds num- 
bers are presented in figure 9. The plots arbitrarily terminate at the body station where 
the two pressure distributions appear to converge. These results indicate that the pri- 
mary effect of displacement thickness on surface pressure is generally confined to the 
first 4 percent of the body length. 

Effect on skin friction .- Skin-friction coefficients for each body with and without 
boundary- layer displacement effects are presented in figure 10 for the maximum and 
minimum test Reynolds numbers. These comparison plots were terminated, in general, 
prior to 1 percent of the body length because the effects were negligible beyond this region. 
This small effect is reflected in the less than 3-percent difference in the integrated values 
of viscous drag with and without displacement- thickness effects, as discussed in the next 
section. The distributions of the skin-friction coefficient, including the effect of boundary- 
layer displacement thickness are shown in figure 11. The location of transition used in 
the calculations is also indicated in figure 11 where appropriate. 

Effect on drag. - The pressure and skin-friction distributions were integrated over 
the entire wetted surface to get the total drag for each body and each test Reynolds number. 
The change in drag coefficient due to boundary-layer displacement effects is observed to 
be nearly constant with power-law displacement at all Reynolds numbers (fig. 12), thereby 
indicating that boundary-layer displacement does not have a strong influence in determin- 
ing the power-law body with the minimum zero-lift drag. The primary effect of displace- 
ment thickness is on the wave-drag coefficient (fig. 12(b)). This effect is more clearly 
evident in figure 13 where the changes in the total-drag-coefficient components due to 
displacement effects are plotted against the power-law exponent for the four Reynolds num- 
bers. The increase in drag due to boundary-layer displacement ranges from 2 to 9 per- 
cent for total drag, 3 to 16 percent for wave drag, and below 3 percent for skin-friction 
drag. The level of the effect is, in general, an inverse function of Reynolds number. 
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Boundary- layer displacement effects on total drag reach a maximum at n near 
0.5 (fig. 13(a)); whereas for wave drag they peak near a value of 0.6 (fig. 13(b)). The 
reduction in displacement effects beyond this peak is more abrupt for wave drag than for 
total drag. In both cases, however, the drop is tempered somewhat by increasing Reynolds 
number. Cross plots of the drag coefficients against Reynolds number (fig. 14) tend to 
obscure these anomalies and, in general, show an inverse variation of boundary-layer 
displacement effects with Reynolds number. The large increases in displacement effects 
for n = 1.00 at the higher Reynolds numbers is attributed to the earlier transition on 
that body. 


CONCLUDING REMARKS 

Wave and skin-friction-drag coefficients have been numerically calculated for a 
series of power-law bodies at a Mach number of 6 and Reynolds numbers, based on body 
length, from 1.5 x 10® to 9.5 X 10®. These calculations were used to evaluate the impact 
of boundary-layer transition and displacement thickness in the determination of the power- 
law body with minimum drag. The calculated-drag results were compared with previously 
obtained experimental data. 

The results of this investigation indicate that, when boundary -layer transition occurs, 
the body for minimum zero-lift drag can be blunter than when the boundary layer is fully 
laminar. This condition occurs because the location of transition at a given Reynolds num- 
ber varies directly with nose bluntness, thereby causing the skin friction to decrease with 
nose bluntness and the minimum total drag to shift to a blunter body. The calculated val- 
ues agree well with the .experimental data. 

The effect of displacement thickness on drag is significant but is nearly constant for 
all the bodies and therefore has very little influence in the determination of the body with 
minimum zero-lift drag. For the present tests, boundary- layer displacement-thickness 
effect is as high as 9 percent of the total drag. 

The results from this investigation are applicable to the design of high-performance 
missiles and cruise and entry configurations at hypersonic speeds. For example, if at 
zero-lift, conditions for boundary-layer transition are encountered in flight, a blunter nose 
should be considered as a trade-off for improving aerodynamic performance as well as 
reducing aerodynamic heating. Furthermore, the boundary-layer displacement thickness 
is a sizable percentage of the body radius, varying from 4 to 12 percent at the base. The 
effect of boundary-layer displacement thickness on drag can also be sizable (up to 9 per- 
cent from this investigation) and should be included in performance calculations. 
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These results are also useful as a guide in the design of hypersonic missiles and 
cruise and entry vehicles for lifting conditions. Previous work has shown that for some 
cases body shapes for minimum wave drag at zero lift have the best performance under 
lifting conditions at hypersonic speeds. Analogously, the identification, in this paper, of 
blunt bodies for minimum drag at zero lift when the boundary layer is transitional should 
be applicable to the design of lifting cruise vehicles and entry configurations. For such 
applications, blunting for practical packaging and aerodynamic stability requirements 
might be less penalizing than indicated by previous studies. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , June 25, 1974. 
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APPENDIX 


DETERMINATION OF BOUNDARY-LAYER DISPLACEMENT EFFECTS 
ON WAVE AND SKIN-FRICTION DRAG 

Procedure 

The effect of boundary-layer displacement thickness on drag was determined by an 
iterative procedure diagramed in figure 15. Seven different computer programs were 
utilized with the iteration requiring three passes through five of the programs to obtain a 
solution. 

To start the calculations, the inverse method for blunt-nose bodies at hypersonic 
speeds (ref. 15) was used to compute the surface pressures on the nose of a sphere and a 
paraboloid and the flow properties in the shock layer along a line normal to each of these 
nose shapes in the transonic flow region. This line is called a field data line in the refer- 
ence. The sphere and paraboloid were fitted to the particular power-law body (n < 1.00) 
at the body slope of the field data line. The flow -field properties of the nose shape, which 
best matched that of the body nose shape, were used in the me tliod- of -characteristics pro- 
gram of reference 16 to compute the pressure distribution over the remainder of the body. 

A curve -fit program using cubic spline functions was then applied to the complete pressure 
distribution. The coefficients of the equation from this pressure curve fit, the body coordi- 
nates from a body geometry program, and the measured location of boundary -layer transi- 
tion were input to the nonsimilar-boundary- layer program of reference 17, which used the 
numerical technique reported in reference 18. In the latter program, the integrated wave 
and skin-friction drags as well as the boundary-layer properties were determined, and the 
coordinates for the body plus the displacement thickness were obtained. 

The iteration procedure for displacement effects began with a program which deter- 
mined an effective power-law body for the basic body enlarged by the boundary -layer dis- 
placement. This effective body was required because the boundary -layer displacement 
was somewhat nonuniform with transition present, and the method-of-characteristics 
program is difficult to use without a continuous surface slope. From this point, the cal- 
culation procedure looped back through the sphere or paraboloid nose match, the method- 
of-characteristics, the pressure curve-fit, and the boundary-layer programs. This loop 
was repeated until the velocity at the edge of the boundary layer agreed with that of the 
previous pass to within 1 percent. Agreement was accomplished in all cases after the 
third pass through the boundary -layer program. 

For the body with n = 1.00 (cone), the initial pass used a constant surface pressure 
from the cone tables (ref. 19), the geometric parameters, and the measured transition 
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APPENDIX - Continued 


location in the boundary -layer program. In the subsequent iterations, the pressure distri- 
butions were obtained from a method-of- characteristics program, developed by Lillian R. 
Boney of NASA Langley Research Center, for the effective body with a conical nose. 

Program Applicability 

Blunt-body pressure. - The only exact nose match for the pressure calculations was 
the paraboloid to the body with n = 0.50. For all bodies with 0.50 g n < 1.00, the parab- 
oloid was the closest match; and the properties for the field data line from it were used 
in the method-of- characteristics calculation. The power-law body with the exact nose 
match to a sphere is n = 0.413; therefore, the spherical nose properties were used only 
for the pressure calculations for n = 0.25. The following table shows the maximum axial 
distances (normalized by body length) between the body nose and the nose of the parab- 
oloid or sphere, which was matched to the body at the initial point of the method-of- 
characteristics solution for the initial pass, first iteration, and the second iteration: 


n 

Ax/Z for - 

Initial pass 

First iteration 

Second iteration 

0.25 

-0.00330900 

-0.00323000 

-0.00323600 

.50 

.00000004 

.00003130 

.00003220 

.60 

.00000841 

.00010470 

.00014150 

.667 

.00003770 

.00005440 

.00005340 

.75 

.00000432 

.00001010 

.00000973 


In all cases, the differences between the actual and calculated nose shapes is 
extremely small. All of the shapes fitted to the power-law bodies lie inside of the body 
except for the sphere fitted to the body with n = 0.25, which is somewhat blunter than a 
sphere. 

Method- of -characteristics solution .- The method-of-characteristics solution extended 
the full body length for all cases except the body with n = 0,75, which covered only 5 per- 
cent of the body. However, comparison of this calculation with the pressure distribution 
obtained with the method-of-characteristics program of Lillian R. Boney, using a 45° start- 
ing cone, showed that the two became coincident at 0.007 percent of the body length from 
the nose; therefore, the two solutions were joined at this point to obtain the complete 
distribution. 

Cubic spline curve fit of pressure distribution .- The cubic spline curve fit was per- 
formed on a cathode ray tube and involved shifting spline curve junction points until the 
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complete curve system was observed to conform to the pressure distribution. In the curve 
fits, from 9 to 16 junction points were utilized, depending on the slenderness of the body. 
Figure 16 shows some examples of the curve fit with the junction points identified by the 
vertical tick marks. 

Boundary-layer program.- The numerical technique presented in reference 18 was 
used to solve the nonsimilar equations governing the laminar, transitional, and turbulent 
compressible flow over the power-law bodies. The turbulent boundary layer was repre- 
sented by a two-layer concept with appropriate eddy-viscosity models used for each layer 
to replace the Reynolds stress term. A static turbulent Prandtl number model was used 
to relate the turbulent-heat- flux term in the energy equation to the Reynolds stress model. 
The mean properties in the transitional boundary layer were calculated by multiplying the 
eddy viscosity by an intermittency function based on the statistical production and growth 
of the turbulent spot in the transitional region of flow. The computer code and flow charts 
are presented in reference 17. 

Effective-body program (Body + 6*).- Before the procedure repeats, the last program 
determines an effective power -law body for the basic body with the boundary -layer dis- 
placement thickness superimposed. Several examples of this curve fit are shown in fig- 
ure 17. The accuracy of representation varies inversely with nose bluntness; nonetheless, 
the representation in all cases is good. 
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Figure 1.- Power-law bodies studied. (All dimensions are in cm.) Model shape equation; 
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Figure 2.- Typical time history of phase change relative 
to local conditions of boundary layer. 
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Figure 4.- Average measured location of initial boundary -layer transition on 
n- powered bodies at Mach number of 6 using phase-change coatings. 
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Figure 5.- Ratio of average measured location of end of transition to initiation 
of transition on n-powered bodies at Mach number of 6. 
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Calculated 



Figure 6.- Variation of calculated skin-friction, wave, and total-drag coefficients with 
power-law exponent in model shape equation — • Measured transition loca- 

tions and boundary- layer displacement-thickness effects are included. 
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(c) n = 0.50; ^ = 1.5 X 10^. 

Figure 9.- Continued. 
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(d) n=0.50; ^ = 9.5 x 10®. 

Figure 9.- Continued. 
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(f) n = 0.60; 9.5 X 106. 


Figure 9.- Continued. 
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(g) n = 0.667; = 1.5 X 10®. 

Figure 9.- Continued. 
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(k) n = 1.00. 
Figure 9.- Concluded. 
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(a) n = 0.25. 

Figure 10.- Effect of boundary-layer displacement thickness on skin-friction 
distribution of each body at minimum and maximum Reynolds numbers. 
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(c) n = 0.60. 
Figure 10.- Continued, 








(e) n = 0.75. 
Figure 10.- Continued 


44 






10 


-1 


10 


-2 





0 .1 .2 .3 A .5 .6 .7 .8 .9 1.0 

S/i 


(a) n = 0,25. 

Figure 11.- Variation of local skin-friction coefficient along body surface for 
four Reynolds numbers. Tick marks locate beginning of transition. 
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(b) n = 0.50. 
Figure 11.- Continued. 
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(a) Total drag. 

Figure 13.- Percent increase in calculated- drag coefficients due' to boundary -layer displacement effects. 
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(a) Total drj 

Figure 14.- Percentage increase in drag coeffi 
layer displacement thickness as fum 






Figure 15.- Calculation procedure for iteration of boundary-layer displacement- 
thickness effects on wave and skin-friction drag for power-law bodies. 























(a) n = 0.25 

Figure 17.- Examples of effective power-law body representations of 
basic body plus boundary-layer displacement thickness. 
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